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ABSTRACT 

Aims. The study of the ensemble properties of the UV/optical broadband variability of quasars is hampered by the combined effects 
of the dependence of variability on timescale, rest-frame wavelength, and luminosity. Here, we present a new approach to analysing 
the dependence of quasar variability on rest-frame wavelengths. 

Methods. We exploited the spectral archive of the Sloan Digital Sky Survey (SDSS) to create a sample of over 9000 quasars in the 
Stripe 82. The quasar catalogue was matched with the Light Motion Curve Catalogue for SDSS Stripe 82 and first-order structure 
functions were computed from the lightcurves. The structure functions are used to create a variability indicator that is related to the 
same intrinsic timescales for all quasars (about 1 to 2 yr in the rest-frame). We study the variability ratios for adjacent SDSS filter 
bands as a function of redshift. A quantitative interpretation of these relations is provided by comparing with the results of simple 
Monte Carlo simulations of variable quasar spectra. 

Results. We confirm the well-known dependence of variability on time-lag; the best power-law fit of the sample-averaged structure 
function has a slope p = 0.31 ± 0.03. We also confirm that anti-correlations exist with luminosity, wavelength, and redshift, where 
the latter can be fully explained as a consequence of the former two dependencies. The variability ratios as a function of redshift 
resemble the corresponding colour index-redshift relations. While variability is almost always stronger in the bluer passband than in 
the redder, the variability ratio depends on whether strong emission lines contribute to either one band or the other. We find that the 
observed variability ratio-redshift relations are described well assuming that (a) the r.m.s. fluctuation of the quasar continuum flux 
follows a power law cr(f A ) oc A" 1 (i.e., is bluer when brighter) and (b) the variability of the emission line flux is only ~ 10% of that 
of the underlying continuum. These results, based upon the photometry of more than 8000 quasars, confirm the previous findings by 
Wilhite and collaborators (2005) for 315 quasars with repeated SDSS spectroscopy. Finally, we find that quasars with unusual spectra 
and weak emission lines tend to have less variability than conventional quasars. This trend is the opposite of that expected from the 
dilution effect of variability due to line emission and may be indicative of high Eddington ratios in these unusual quasars. 

Key words. Galaxies: active - Galaxies: statistics - Quasars: general - Quasars: emission lines 



1. Introduction 

The variation in flux density over time is one of the key charac- 
teristics of active galactic nuclei (AGN). From the onset of the 
study of quasars, which are the most luminous AGNs, variability 
is known to be a diagnostic property and has been successfully 
used as a criterion for selecting quasar candidates in a number 
of studies. On the other hand, the physical mechanisms behind 
these fluctuations are still poorly known, though it has long been 
understood that the observed flux variations of AGNs hold keys 
to the structure of the radiation source. The advent of modern 
variability surveys (e.g., Ivezic et al. 2003; Groot et al. 2003; 
Becker et al. |2004l Bauer et al. |2009l Djorgovski et al. |2008l 
LSST Science Collaborations 2009) combining high photomet- 
ric accuracy, large survey areas, sufficiently long time baselines, 
and a high number of observation epochs holds the promise of 
significant progress. 

Numerous studies have shown that the optical/UV contin- 
uum variability of quasars correlates with various physical pa- 
rameters such as time-lag, rest-frame wavelength, luminosity, ra- 
dio properties, emission line properties, and the Eddington ratio 
(e.g., UMchet al. [T997l Giveon [T999l H elfand[200D Gas kell & 
Klimek|2003j Vanden Berk et al. 120041 Rengstorf et al. 120061 



W old et al. 120071 Wilhit e et al . 120081 Bauer et al. [20091 Ai et 
al. 120101 Schmidt et al. 1201 01) . Despite considerable progress 
in the description of quasar variability, there are still conflict- 
ing claims about these correlations and their physical interpre- 
tation is fraught with problems. The most frequently discussed 
models for the continuum variability of AGNs include instabil- 
ities in the accretion disk (e.g., Kawaguchi et al. [1998), mul- 
tiple supernovae (or other Poissonian processes; e.g., Terlevich 
et al. |1992t Cid Fernandes et al. 119971 1, star-star collisions in 
the dense circumnuclear environment (Torricelli-Ciamponi et al. 
2000), and microlensing by compact foreground objects (e.g., 
Chang & Refsdahl [T979l Ha wkins [T993H20T01 Lewis & Irwin 
[T9961 Zackrisson et al. 120031 . 

Much of the knowledge about quasar variability is based 
on the detailed investigation of large, statistically well-defined 
quasar samples. Even though each individual lightcurve is 
poorly sampled, a large number of quasars and the wide range 
of parameter values allow us to distinguish the dependences 
of variability on the various parameters. Hence, substantial im- 
provement has been achieved by employing the quasar data base 
from the Sloan Digit al Sky Survey (SDSS; York et al. [2000; 
Abazajan et al. 120091 1. Using a sample of 25 000 spectroscop- 
ically confirmed quasars from the SDSS, Vanden Berk et al. 
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(2004) parametrized correlations between variability and time- 
lag, luminosity, redshift, and rest-frame wavelength upon a time- 
baseline of up to ~ 2yr in the observer frame. A sample of 
more than 41 000 quasars was studied by de Vries et al. ((2005 ). 
By combining the SDSS photometric data with the Palomar 
Observatory Sky Survey (POSS), these authors achieved a base- 
line of up to 50 yr in the observer frame, yet with a very sparse 
lightcurve sampling per quasar. 

Wilhite et al. (2005 ) selected a set of 315 quasars with multi- 
epoch SDSS spectroscopy to analyse variability from spec- 
trophotometry. In this way, the wavelength dependence of quasar 
variability could be studied at spectral resolution high enough 
to allow the analysis of variability not only in the continuum 
but also in the major emission lines. The results confirm that 
the continua are bluer when brighter and clearly show that vari- 
ability is weak in the emission lines. They provide a natural ex- 
planation of the intrinsic Baldwin effect behaviour of individual 
AGNs (Kinney et al. [T9901 Wilhite et al. I2006I Pereyra et al. 
(2006) demonstrated that the composite flux differences in the 
rest-frame wavelength range 1300-6000 A can be fit by a stan- 
dard thermal accretion disk model in which the accretion rate 
has changed from epoch to epoch. These results suggest that 
the dominant fraction of the optical/UV variability of quasars 
is related to the accretion process. On the other hand, chro- 
matic variability may also be produced by quasar microlensing 
as the continuum source is compact and bluer at smaller radii 
(Wambsganss & Paczvhski ll9911 Yonehara et al. 2008 and ref- 
erences therein). Wilhite et al. ((2005 ) argued that the low vari- 
ability of the broad emission lines provides a problem for the 
microlensing scenario as one would expect a higher fraction of 
events where the lense crosses the broad line region (BLR) but 
not the compact source. The BLR was once considered too large 
to be affected significantly by microlensing. Studies have shown, 
however, that the BLR can be small enough to be amplified sig- 
nificantly by stellar-mass objects (Abajas et al. 2007, and refer- 
ences therein). 

A 300 square degree area in the Southern Galactic Cap has 
been repeatedly imaged by the SDSS since 1998 and by the 
SDSS-II Supernova Survey (Frieman et al. 120081 ) since 2005. 
The multi-epoch observations in this SDSS Stripe 82 (S82) pro- 
vide a useful database for the statistical analysis of quasar vari- 
ability. Sesar et al. (2007 ) were the first to analyse the variable 
sources in S82 with the emphasis of characterizing the faint vari- 
able sky and quantifying the variable population as a whole. 
They present methods to select variable sources and discuss their 
distribution in the magnitude-colour- variability space. The three 
dominant classes of variables found by Sesar et al. are quasars 
(63%), RR Lyrae stars (7%), and stars from the main stellar lo- 
cus (25%). Moreover, they demonstrate that at least 90% of the 
quasars are significantly variable. An interesting detail of this 
paper is a plot (their figure 7) of the ratio, V g /V r , of the variabil- 
ity V of quasars in the g- and r-band as a function of redshift z. 
Sesar et al. observe a feature in this diagram at z ~ 1.0-1.6 and 
argue qualitatively that it is due different parts of the spectrum 
varying in different ways. This is an important observation that 
offers a new possibility of investigating the spectral variability 
of quasars, which is obviously worth studying in detail. 

The present paper is focussed on the variability ratios 2k,i = 
Vk/Vi as a function of redshift. The main aim is to derive the 
dependence of the variability on the (intrinsic) wavelength. In 
view of the combined effects of the dependence of the vari- 
ability on wavelength, time, and luminosity, the ratios of the 
variabilities in different photometric bands provide a better ap- 
proach to the spectral behaviour of variability than variability 



measured in a single band, as was pointed out by Sesar et al. 
(120071 ) . The main differences from Sesar et al. are the following: 
(i) Sesar et al. used photometric measurements from 58 imag- 
ing runs from 1998 September to 2004 December. In the present 
study, the lightcurves are taken from the Light and Motion Curve 
Catalogue (LMCC) of the SDSS S82 (Bramich et al. l2008l . The 
catalogue is based on 134 imaging runs obtained from 1998 to 
2005 November with 62 runs obtained in 2005 alone, (ii) We ex- 
tend the approach to all five SDSS bands, (iii) We use variabil- 
ity indicators related to (more or less) the same quasar-intrinsic 
timescale to eliminate z-dependent selection effects due to cos- 
mological time dilation, (iv) Since the approach works only if 
quasar variability is dominated by the same processes at low 
and high redshift, we check whether this assumption is consis- 
tent with our data, (v) Finally, a quantitative interpretation of 
the diagrams is presented by means of numerical simulations. 
Compared to the spectroscopic sample of Wilhite et al. ((2005 ), 
the quasar sample of the present study is much larger and the 
time baseline covered by the data is longer. 

The creation of the quasar sample is described in Sect. [2] 
Variability indicators are computed from their lightcurves in 
Sect. [3] The resulting Q^-z relations for the five SDSS bands 
(k,l — 1 ... 5) are presented in Sect.|4]and compared to the re- 
sults of a simple parametrized model in Sect. [5] Other subjects 
in Sect.[5]are the comparison of our results with Wilhite et al. 
(2005 ) followed by a brief discussion of the variability proper- 
ties of unusual quasars and the subsample of radio-loud quasars. 
Finally, Sect.[6]presents our conclusions. Standard cosmological 
parameters Hq = 71 km s _1 Mpc~',Q m = 0.27, Qa = 0.73 are 
used throughout the paper. 

2. The quasar sample 

The SDSS Seventh Data Release (DR7; Abazajan et al. 120091 , 
marking the completion of the original goals of the SDSS and 
the end of the phase known as SDSS-II, contains over 1.6 mil- 
lion spectra in total, including 120 000 quasars. The quasar sam- 
ple for the present study was constructed from the SDSS DR7 
public archive of spectra in the region of the Stripe 82, i.e., right 
ascension a = 20 h . . .4 h and declination 6 = -1°25 . . . + T25. 

Whether a spectrum is useful for scientific purposes or not 
is indicated by the SciencePrimary flag, which is set to be ei- 
ther 1 or 0. Since SciencePrimary = 1 does not necessarily 
guarantee the correctness and precision of the cataloged quan- 
tities, we decided to download all spectra classified as quasars 
(Spec_cln=3), regardless of their SciencePrimary flag, including 
all spectra classified as unknown (Spec_cln=0). This last deci- 
sion was motivated (i) by our aim of creating a quasar sample as 
large and complete as possible and (ii) the claim that interesting 
cases of rare unusual quasar spectra (e.g., Hall et al. 2002 ) might 
be hidden there. The list contains 20 542 objects of which 19 660 
were available for download (13 224 with Spec_cln=3 and 6436 
with Spec_cln=0). The quasar selection does not include AGNs 
classified spectroscopically as galaxies (Spec_cln=2). We expect 
the majority of them to be low-redshift, low-luminosity AGNs, 
while the present study aims to compare the variability at dif- 
ferent redshifts z = 0.2 ... 3 and is restricted, therefore, to high- 
luminosity AGNs (M; < -21, Sect.|4.3|). Moreover, the stellar 
light from the host galaxies of low-luminosity AGNs tends to 
"dilute" the variability of the AGN. 

All spectra were checked individually to estimate redshifts 
z and object types. Redshift values were computed automati- 
cally by comparing of measured line centres of the prominent 
emission lines with a catalogue of their rest-frame wavelengths. 
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The spectral lines useful for the redshift (z) estimation were se- 
lected manually: this allows us to recognize and exclude quasar 
lines influenced by e.g., the remnants of bright night sky lines 
as well as Lya lines affected by strong absorption of the blue 
wing. Dubious cases were flagged (Flag_z = 1). While the ap- 
proach was not optimized to provide the highest accuracy z val- 
ues, roughly erroneous redshift estimations could be efficiently 
suppressed, which is more important for the present study. The 
spectra were classified into the six types: quasar (q), possible 
quasar (q?), galaxy (g), star (s), unclear (u), and noise (n). For 
quasar spectra, special features such as strong broad absorption 
lines (BALs), unusual BALs, weak emission lines, or strong iron 
emission were noted. The procedure was performed twice for all 
19660 spectra. 

Table 1. Number of quasars from the basic catalogue with red- 
shift deviations |z here - zsdssI > |Az| max . 



|AzUx= 0.05 0.1 0.5 1.0 2.0 




1 2 3 4 5 lOOO 

rsdsriift z t per run 



Fig. 1. Left: Redshift versus absolute i band magnitude and his- 
tograms for the distributions of redshifts {top) and absolute mag- 
nitudes {right) for the final quasar sample. 



8 212 quasars found in the QC DR5 : 

N{> |AzU) 22 5 3 1 

1 643 quasars not found in the QC DR5: 

N{> ]Az| max ) 391 381 314 167 51 



We assigned 13 558 selected spectra to the types q, q?, or g. 
To distinguish galaxies from quasars, a simple criterion was used 
where priority was given to low contamination by non-AGNs: 
spectra containing typical absorption lines and/or narrow emis- 
sion lines but no signs of broad emission-line components were 
classified simply as galaxies. The total number of quasar spec- 
tra (type q) is 12 966, among them 976 from the database of 
the unknowns (Spec_cln=0), another 4 objects were classified as 
possible quasars. Fiber positions that differ by less than the fibre 
diameter (3 arcsec) can be ascribed with high probability to the 
same object. After excluding multiple observations, the resulting 
basic catalogue contains 9 855 entries with 2021 objects having 
been observed at least two times (660 objects with more than 2 
spectra). For all objects with multiple spectra, the z and type of 
the single spectra were compared with each other. Inconsistent 
results were found in only very few cases (e.g., because of cor- 
rupted spectra or low signal-to-noise ratio (S/N)). 

We identified 8 212 entries in our basic quasar catalogue with 
entries in the SDSS DR5 quasar catalogue (QCDR5; Schneider 
et al. 12007b . 736 (7%) were found only in the database of the un- 
knowns only. Table[T]compares our redshift estimates with the z 
values from Schneider et al. (2007 ) and the DR7, respectively, 
for objects not listed in the QCDR5. The 25 quasars with dubi- 
ous redshift flags from the present study were eliminated. While 
the agreement with the QCDR5 is satisfactory, we find that a 
substantial number of erroneous redshifts have been produced 
by the SDSS DR7 pipeline, which are attributed to erroneous 
line identifications (e.g., Mg II with Lya, C IV, or C III]) in most 
cases. 

To analyse variability, the basic quasar catalogue was 
matched to the LMCC (Bramich et al. 2008) using an iden- 
tification radius of 2 arcsec. The LMCC was trimmed to only 
those objects that have a mean right ascension in the range 
a = 20. 7 h . . . 3.3 h because of the sparse temporal coverage out- 



side these limits. Compared to the whole S82, the area of sky 
covered by the LMCC is thus reduced by ~ 17% to ~ 249 deg 2 . 
The final quasar catalog contains 8 744 objects identified in the 
LMCC (Fig.[T}. The redshifts cover the range z ~ 0.1 ... 5 with 
Z = 1.46 + 0.79. Absolute magnitudes M; were computed fol- 
lowing Kennefick and Bursick (2008; their equations 1 to 6) 
correcting the SDSS i band magnitude for Galactic extinction 
and applying a standard continuum K correction by assuming 
a spectral index a v - -0.5. The Galactic foreground extinction 
was obtained from the "Galactic Dust Extinction Service'Q of 
the NASA/IPAC Infrared Science Archive using the translation 
relations for the SDSS bands given by Schlegel et al. ( 119981 ). 
The mean value of the absolute magnitudes in the sample is 
M; = -24.7 + 1.57. The distribution of the quasars in the z-M\ 
plane is shown in Fig.[T] 

To model spectral variability (Sect. |5.lj ), we needed to pro- 
duce an average quasar spectrum. We used our "standard sam- 
ple" (see Sect. 3.1 below) and constructed a composite follow- 
ing Vanden Berk et al. (2001). As we are primarily interested 
in "normal" quasars, objects with highly unusual spectra (re- 
mark = 'x_bal' or 'myst' in Table|2]) are not included as well 
as quasars with uncertain redshifts. We also note that the sam- 
ple is corrected for multiple observations; quasars with two or 
more SDSS spectra are thus represented by only one spectrum, 
usually the one with the highest S/N. All spectra were resized to 
a common range of 3900-9 100A in the observer frame, sorted 
by increasing redshift, corrected for Galactic extinction, shifted 
into their rest-frame, and binned to a resolution of 1 A per pixel. 
Subsequently each spectrum was inserted as a new row into a 
two-dimensional image, with 1A binning on the horizontal axis 
and one spectrum per pixel on the vertical axis. Every time be- 
fore a new spectrum was inserted, a preliminary mean quasar 
spectrum was created simply by computing the mean pixel value 
for each wavelength bin in the two-dimensional image. The re- 
sult is used for the normalisation: before inserting into the two- 
dimensional image, each spectrum is normalised to the mean 
pixel value of the overlapping region with the respective prelim- 
inary composite spectrum. The spectra was not assigned any sta- 
tistical weights, i.e. the S/N of the spectra is assumed to be rela- 
tively constant throughout the sample and over the observed win- 



http://irsa.ipac.caltech.edu/applications/DUST/ 
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Fig. 2. Top: Median composite spectrum from the final quasar 
sample with fitting power-law continuum. Second row: The 
semi-interquartile range, SIQR, in units of the median. Third 
row: Relative deviation from the Vanden Berk et al. (2001) com- 
posite ([here-Vanden Berk]/here). Bottom: Number of spectra 
per wavelength bin used for the construction of the composite. 



dow. Pixels flagged 'bad' in the SDSS 2D spectroscopic pipeline 
are interpolated over, and outliers and negative fluxes were ex- 
cluded. 

Composite spectra can be created in different ways. Vanden 
Berk et al. (2001 ) discussed how arithmetic and geometric mean 
composites have different properties. Following the arguments 
given by Wilhite et al. (2005 ), we consider the arithmetic com- 
posite but we prefer to use the median instead of the mean. The 
final composite spectrum is constructed by computing the me- 
dian value for all wavelength bins represented by at least three 
spectra. The result is shown in the top panel of Fig. [2] 

Vanden Berk et al. (JJOOl) suggested using the semi- 
interquartile range, SIQR= (q^ - q\)l2 as a measure of the 
spectrum-to-spectrum variation, where q\ and qj are the first 
(25%) and third (75%) interquartile; the second interquartile, qi, 
is the median. This ratio, SIQR/<?2, is shown in the second panel 
of Fig. [2] There exist many spectra with a relative flux much 
higher or lower than the median at wavelengths where SIQR/<72 
is high . On the other hand, if SIQR/^2 is small, the distribu- 
tion is sharply peaked. As expected, the variation is huge in the 
wavelength range of the Lya forest. At longer wavelengths, the 
variations are strong in the peaks of the emission lines, in par- 
ticular of the [Om] lines. Longward of ~ 4000A the spectrum- 
to-spectrum variation increases because of the different strong 
contribution of the host galaxies. 



The third panel of Fig. [2] compares our median compos- 
ite with that of Vanden Berk et al. (2001); where both spectra 
have been normalized normalized at 4000 A. For wavelengths 
A - 1000 . . . 4000 A, our result (Fig.||]) is in accord with the 
Vanden Berk composite. The deviation at longer wavelengths re- 
flects the differences in the quasar selection for the two samples. 
Vanden Berk et al. define "quasar to mean any extragalactic ob- 
ject with at least one broad emission line and that is dominated 
by a nonstellar continuum". Spectra with continua dominated by 
stellar features were rejected by these authors "to avoid intro- 
ducing a significant spectral component from the host galaxies". 
This criterion was not applied in the present study. The con- 
sequence is a stronger contribution from host galaxies in our 
spectrum, yet the absolute magnitude-redshift plot appears to 
have a similar representation from low-luminosity AGNs in both 
cases. At short wavelengths, the differences are likely statistical 
fluctuations due to the small number of spectra included (bot- 
tom panel) in combination with the strong spectrum-to-spectrum 
variations at these wavelengths. In the spectral region between 
Lya and 4000 A, the continuum is matched by a power law 
fx cc A aA with ax = -1.52 (a v = -0.48, in good agreement 
with a v = -0.46 given by Vanden Berk et al. for their arithmetic 
median composite). 

The quasar catalog was matched to the FIRST (Becker et 
al. 1 19951) VLA 20 cm catalogue (match radius 2 arcsec). For the 
467 identified FIRST sources, the values for the integrated 20 cm 
flux and the peak flux were entered. A radio compactness pa- 
rameter and the radio loudness parameter R[ = log (,F ra dio/f i) 
were computed following Ivezic et al. (2002); 381 sources are 
radio-loud (/?; > 1). Finally, more than 50 quasars were marked 
as having more or less unusual spectra, i.e., spectra that devi- 
ate remarkably from the typical quasar spectrum. The classifica- 
tion was performed by eye and is thus not based on quantifiable 
properties. The spectra of these objects display various unusual 
features that are difficult to discern at first glance, often (but 
not exclusively) related to complex broad absorption line (BAL) 
troughs combined with strong reddening (e.g., Hall et al.2002). 
This subsample is inhomogeneous and includes unusual BAL 
quasars as well as weak emission line quasars (WLQs) and 
"mysterious" objects such as SDSS J010549.7 5-0033 14.0 and 
SDSS J220445.26+003 142.0 from Hall et al. ( [20021 . In sev- 
eral cases, the redshifts are uncertain. For one of these objects, 
SDSS J014349. 14+002128.4, the spectrum was first classified as 
an unusual BAL quasar but later rejected after we found it to be 
a superimposition of the spectra of a normal quasar in the blue 
and an M star in the red. 

3. Measuring spectral variability 

3.1. Photometric data 

The catalogue of objects in S82 compiled by Bramich et al. 
(2008) contains about 3.7 million stellar objects and galaxies. 
It comes in two flavours, the Light-Motion Curve Catalogue 
(LMCC), which contains the measured quantities for each ob- 
ject listed as a function of filter band and epoch, and the Higher- 
Level Catalogue (HLC), which contains similar quantities for 
each light-motion curve. Since the present paper intends to ana- 
lyze the lightcurves, the LMCC is used here. For the overwhelm- 
ing majority (<: 90 %) of the quasars, the number of measure- 
ment epochs is between 20 and 50 with a mean at ~ 35 in all 
five bands. For a detailed description of the LMCC, we refer to 
Bramich et al. (2008 ) and address only three issues: (i) for an 
object record to be included in the LMCC, a list of constraints 
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need to be satisfied in at least one waveband. In addition, we ac- 
cept here only measurements, in each band, with the quality flags 
BRIGHT=0, EDGE=0, BLENDED=0, and SATUR=0. (ii) As 
noted by Bramich et al., the catalogue is contaminated by some 
strong photometric outliers that are clearly due to 'bad epochs'. 
We find that some lightcurves show one, sometimes two, rarely 
more, examples of such low quality data that can significantly 
affect the measurement of variability. We applied a 2.5<xioo- 
clip criterion to identify and eliminate the corresponding data 
points, where o"ioo is the standard deviation of magnitudes in 
a 100 d time interval centred on the corresponding epoch. Each 
lightcurve was checked for up to three of these outliers, (iii) We 
restrict the analysis, in each band, to objects brighter than the 
limiting magnitudes 21.9,23.2,22.8,22.2, and 20.7 for u,g,r,i, and 
Z. More precisely, only quasars with mean magnitudes 0.5 mag 
brighter than the corresponding limit in each band were consid- 
ered to minimize a possible bias by a one-sided censorship of the 
lightcurves due to variability. 

As a first step, we estimate the fraction of significantly vari- 
able quasars. We use the x\- a P er degree of freedom defined by 
Sesar et al. (2007 ) to compare the magnitude fluctuations of the 
quasars from the final sample with those of 1 100 objects se- 
lected from the SDSS standard star catalogue for S82 (Iveciz et 
al. 120071) . The constraint x\ qqs > ^ (significance level a = 0.05) 
is a useful criterion for significant quasar variability, i.e., the 
fluctuations cannot be explained by Gaussian-distributed mea- 
surement errors with the same (magnitude-dependent) standard 
deviations as measured for the standard stars. Using this crite- 
rion, we find that 93%, 97%, 93%, 87%, and 37% of the quasars 
are variable in the u,g,r,i, and z band, respectively. 

We next consider whether variability is characterized more 
accurately by a normal distribution of fluxes or magnitudes. 
Fig. [3] shows the distribution of the normalized fluctuations 
Axjj = (xij-Xj)/cr x j for the g-band measurements i of the signif- 
icantly variable quasars j, where x represents either magnitudes 
(left hand side) or fluxes (right hand side) after the exclusion of 
outliers. The data are obviously more accurately matched in the 
left panel, yet the agreement is not perfect, found for both the 
optical and X-ray radiation of AGNs (Gaskell & Klimek I20031 
and references therein). 




2 2 

(if [g]-f [g] )/o f[g] 



Fig. 3. Distribution of the normalized fluctuations of g-band 
magnitudes (left) and fluxes (right), respectively. Dotted curves: 
normalized, centred Gaussian distribution. 



As we are interested in the variability properties of typical, 
'normal" quasars, the unusual BAL quasars and the quasars with 



uncertain redshifts were removed from the sample. We also re- 
jected the core-dominated radio-loud quasars. Vanden Berk et 
al. (|2004) and Rengstorf et al. (|2006) found the median inte- 
grated radio-flux density to be higher for the variable population 
of their samples than for the nonvariable population and sug- 
gested blazars as a possible explanation. The continuum emis- 
sion mechanism for blazars differs fundamentally from that of 
other quasars (e.g., Ulrich et al. 1997). It is unlikely that all re- 
jected radio-loud quasars belong to the blazar family but there 
might be a substantial fraction where both the flux and its vari- 
ation are dominated by the relativistic outflow. In the following, 
we consider as our "standard sample" all quasars with known 
redshifts from the final sample without the unusual and the 
known radio-loud quasars. The latter two subsamples are dis- 



cussed separately in Sect. 5.4 



3.2. Structure functions 

The statistical measure of quasar variability used in the present 
paper is based on the first-order structure function (SF) D(t), 
where r is the time-lag between two observations. The SF is 
closely related to the autocorrelation function and contains, as 
a sort of running variance (as a function of the time-lag), in- 
formation about the timescales of the involved variability pro- 
cesses. The SF was first introduced to quasar research in radio 
astronomy (Simonetti et al. 119851 Hjellming and Narayan 1986; 
Hughes et al. 1992 ) and has also become a popular tool in optical 
studies of both quasar samples and individual quasar lightcurves 
(e.g., Giallongo et al. 119911 Hook et al. 119941 Meusinger et al. 
[19941 Kawaguchi etal. [19981 Vanden Berkl2004l De Vriesetal. 
120051 Wilhite et al. l20Q.Sl Bauer et al. l2009l Schmidt et al. l20T0l . 
Scholz et al. (Q997 ) were the first to use an index for long-term 
variability based on the observed SF to identify quasar candi- 
dates in an variability-based quasar survey (see also Meusinger 
et al. [2002b . 

There are three different definitions of the SF used in the 
literature. Here, we follow the version introduced in the pioneer- 
ing papers of Simonetti et al. (1985) and Hjellming and Narajan 
( 1986), where the first-order SF is defined as 



D(t) = ([s(t + t) — s(t)] ) t , 



(1) 



where s(t) is the measured signal at time t, and the angular 
brackets denote the time-average. The SF analysis in the opti- 
cal/UV usually identify signals with apparent magnitudes, i.e., 
s(t) = m(t). In a second approach (e.g., DiClemente et al. 1996; 
Vanden Berk et al . 120041 Wilhite et al. l2008l ([Am] 2 ) is replaced 
by 7r/2(|Am|) 2 . The former was referred to as SF (A) and the lat- 
ter as SF (B) by Bauer et al. (2009 ). In a third version, here called 
SF (C) , (|Am|) was used instead of ([Am] 2 ) where () means either 
the median (Hook et al. |1994| l or the mean (Schmidt et al. 1201 Oi l. 
Sp(A) nas (jjg ac i van tage of being directly related to other statisti- 
cal quantities such as the variance and the autocorrelation func- 
tion, while SF (B) and SF (C) are less sensitive to the presence of 
outliers in the data. Since we rejected the outliers in the process 
of constructing the lightcurves (Sect. |3.lj ), we decided to apply 
the "classical" definition from Eq.[T] 

The general definition and some properties of the SF are 
given by Simonetti et al. 0985), an "ideal" SF of a "typical" 
measured process is schematically discussed by Hughes et al. 
(II 992l l. For lags shorter than the smallest correlation timescale, 
To, of the variable process, the SF is given by Do = D(j <K To) = 
2£ 2 , where £ 2 is the variance of the measurement noise. The 
SF corrected for the measurement noise is given by D coit (t) = 
D(t) — 2£ 2 . For t » T\, the SF displays another plateau with 
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Table 2. First five rows of the table containing the variability indicators Vt. The full table is available at the CDS. 



number 


RA 


DEC 


redshift 


flag 


g 


Mi 


V u 






Vi 


V z remark 




(degrees) 


(degrees) 










(mag 2 ) 


(mag 2 ) 


(mag 2 ) 


(mag 2 ) 


(mag 2 ) 


1 


0.001997 


-0.451102 


0.250 





21.257 


-20.39 


0.000 


0.034 


0.012 


0.013 


0.004 


2 


0.027231 


0.515332 


1.823 





20.212 


-25.20 


0.073 


0.122 


0.088 


0.027 


0.009 


3 


0.033946 


0.276292 


1.837 





20.213 


-25.75 


0.012 


0.025 


0.017 


0.005 


0.013 


4 


0.049839 


0.040359 


0.479 





17.859 


-24.15 


0.044 


0.017 


0.005 


0.003 


0.000 


5 


0.051084 


-0.539051 


1.436 





20.435 


-24.56 


0.179 


0.165 


0.091 


0.057 


0.080 



D(t) = 2cr 2 , where T\ is the longest correlation timescale and 
cr 2 > £ 2 is the variance of the process. For the intermediate part, 
To ^ t ^ T\, the SF of a stationary random process is character- 
ized by a power law D(j) oc t a , where the exponent depends on 
the power spectrum of the process. 

We briefly note two other fundamental properties of the 
SF. First, the method is suited to the analysis of sparsely sam- 
pled lightcurves of non-periodic and non-sinusoidal processes. 
Second, provided that variability is quasar-intrinsic, the data ob- 
tained at the same epochs for quasars at different z cover dif- 
ferent source-intrinsic timescales because of the cosmological 
time-dilation and are thus not directly comparable. This prob- 
lem can be avoided using the SF as a function of the rest-frame 
time-lag t, = r /(l + z). Giallongo et al. ( 1199H were the first to 
introduce an indicator for optical quasar variability based on the 
rest-frame SF. 

In the following, variability is measured by the first-order 
structure function, corrected for measurement errors and binned 
into rest-frame time-lag intervals 



V = (D( Tr )) Ti - 2?m, 



(2) 



where % 2 (m) is the variance in the measurement noise for star- 
like objects of the apparent magnitude m = m adopted from the 
standard stars, where in is the mean magnitude of the quasar. 
The angular brackets denote the arithmetic mean over a suitable 
interval of rest-frame time-lags around r, . 

3.3. Spectral variability 

For the investigation of the wavelength dependence, the variabil- 
ity is measured by one single value V\ for each quasar in each 
of the five SDSS bands k = 1 ... 5. Distinguishing the V(A, L, z) 
relations requires two steps. As suggested by Sesar et al. (2007 ), 
we consider the ratio Qk.i = Vjc/Vi of the variabilities in two 
adjacent photometric bands (k = 1 . . .4, 1 = 2 ... 5) instead of 
Vk itself. Second, to eliminate z-dependent selection effects due 
to cosmological time dilation, the r r interval has to be chosen 
properly. The latter is constrained by the following criteria: 

- Equation[TJ yields discrete data points for each quasar in the 
D-t diagram. The mean SF {D(t t )) Ti is computed by averag- 
ing these data within discrete r r bins. The binning intervals 
must be wide enough to cover a large number of data points. 

- D(t t ) must be related to the same T t interval for all quasars 
independent of z. 

- Because quasar variability increases monotonically with t> 
up to timescales of ~ 40 yr (de Vries et al. 2005), r r should 
be as long as possible. 

- The maximum rest-frame time-lag is T r max ~ 7 yr (for z ~ 0), 
set by the time baseline of the LMCC. 

- To investigate the wavelength dependence, variability has to 
be measured in narrow z bins (see below). The z distribution 



for our final quasar sample (Fig.[T| shows a strong decline 
at z ^ 2. Beyond z ~ 3, the number of quasars per z bin 
becomes too small for a statistical analysis. 

Taken together, we find that r r = 300 - 600 days is an appropri- 
ate binning interval. We computed (D(r r )) Tr to be the arithmetic 
mean of all D(t v ) in this r r interval. The value Vk( T r) from the 
SF binned in this way is considered a useful variability indicator. 
The present analysis is thus related to variability processes with 
typical timescales of the order of ~ 1 ... 2 years (rest-frame). 
The median rest-frame time-lag of 430 d is a factor of 3.5 longer 
than for the spectroscopic sample of Wilhite et al. (f2005). This 
is an important advantage of the present data because of the in- 
crease in variability with time-lag (Fig.Q. 

The variability indicators Vk for the 8 744 quasars from our 
final sample are given in Table|2] available at the CDS. The 
table contains the following information. Column 1 is the run- 
nung number, Columns 2 and 3 list right ascension and declina- 
tion (J 2000) taken from the fits headers of the SDSS spectra. 
Columns 4 and 5 give the redshift and the redshift flag (=1 for 
uncertain redshifts, otherwise). The g band magnitude is given 
in Col. 6 and the absolute magnitude M\ is given in Col. 7. The 
following five columns contain the variabilities V for ugriz. For 
(D(t)) < 2£ 2 , the value of V\ was set to zero, because negative 
values of Vk are unphysical. The last column contains informa- 
tion on specific spectral features (w_l: weak-line quasar, s_bal: 
strong BAL structure, x_bal: extreme BAL structure, myst: mys- 
terious spectrum). 



4. Results 

4.1. Variability as function of time-lag 

The binned, "noiseless" first-order structure functions averaged 
over our standard sample of quasars with M[ < -22 for u,g,r,i, 
respectively, are shown in Fig.|4] Seven time-lag bins were cho- 
sen in such a way that we have, on the one hand, log to intervals 
of comparable widths and, on the other, reasonably large num- 
bers of measurements in each bin. The vertical bars indicate the 
shift due to the corrections for instrumental errors ,which are 
large on short time-lags. As is evident from Fig.|4] variability in- 
creases with time-lag and decreases with rest-frame wavelength. 
The detailed investigation of the dependence of variability on 
rest-frame wavelength is the subject of Sect. 4.3 Here we dis- 
cuss the dependence on time-lag. 



In each of the four bands, the SF is closely fitted by a power 



(following the notation of Kawaguchi et 



la w [V(y r )] 1/2 oc 
al. U998b . No significant flattening of the SF is observed at long 
time-lags. The slope of V(t) in the double-logarithmic presen- 
tation depends slightly on the chosen lag intervals and also on 
the absolute magnitude threshold for the quasar sample. Results 
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Table 4. SF exponent ft from previous work and the present 
study. 



log <time-lag (days!) 

Fig. 4. Corrected ensemble-averaged structure functions for the 
quasars with M\ < —22 in the u,g,r, and i band as a function 
of time-lag. The time-lags are in the rest-frame, yet the SDSS 
bands are in the observer frame. 

Table 3. SF exponent ft for different parameter selections. 



Reference 




ft 

H 


r mfl , (vr) 


SF 


Hooketal. J 1994 





.20 


6(rf) 


C 


Kawaguchi et al. 11998) 





.35 


2 (of) 


A 


Hawkins 12002&, Quasars 


0.20 


±0.01 


20 (of) 


A 


Hawkins (12002). Seyferts 


0.36 


±0.02 


20 (of) 


A 


Vanden Berk et al. 12004) 


0.25 


±0.01 


2(rf) 


B 


deVries et al. 12005) 


0.30 


±0.01 


50 (of) 


A 


Rengstorf et al. ( 2006 1 


0.47 


±0.07 


2(rf) 


B 


Wilhite et al. 12008) 


0.46 


±0.03 


2(rf) 


B 


Bauer et al. 12009) 


0.36 


±0.01 


3(rf) 


A 


Schmidt et al. 12010) 


0.43: 


t 0.002 


10 (of) 


C 


present study 


0.31 


±0.03 


3(rf) 


A 



M; < -21 



M; < -22 



M; < -23 



Ty > 10 days u 



0.230±0.011 
0.314±0.031 
0.322±0.024 
0.305±0.010 



0.231±0.010 
0.326±0.038 
0.326±0.034 
0.339±0.013 



0.231+0.01 1 
0.331±0.046 
0.368±0.048 
0.373±0.025 



r r > 100 days u 



0.237= 
0.278= 
0.299= 
0.311= 



0.010 
0.010 
0.017 
0.009 



0.236±0.009 
0.279±0.010 
0.312±0.019 
0.331±0.012 



0.238±0.011 
0.272±0.009 
0.322±0.021 
0.345±0.013 



for the value of the exponent ft are listed in Table[3] The un- 
certainties given there are the rms errors in the regression only. 
The total errors are larger and are expected to be dominated by 
the uncertainties in the correction for measurement errors, which 
may be underestimated here especially in the u band. (An under- 
estimation lowers the SF slope.) For M; < -22, the mean value 
of ft, averaged over the high-throughput SDSS g,r, and i filter 
bands, is ft = 0.33 + 0.01 for r r > 10 days and 0.31 ± 0.03 for 
r r > 100 days, respectively. The latter value is probably more 
representative because the SF is less sensitive to the correction 
for measurement errors at longer lags. 

A summary of SF slopes found in the literature is given in 
Table|4] The values for ft scatter far more than expected from the 
individual errors. However, one has to take into account that the 
results are based not only on different quasar samples but also 
on different definitions of the SF (last column; see Sect. |3.2| i, 
different methods of averaging, different corrections of mea- 
surement errors, and different time baselines. Therefore, Table|4] 
also lists the definition and the maximum time-lag r max (where 
"of" and "if" refer to observer frame and rest-frame, respec- 
tively). The simple average of the previous results based on the 
SF (A) , except for the Seyfert value from Hawkins (2002), yields 



0.30 + 0.07 (0.34 + 0.1 1 for all data) which is close to the value 
from the present study. Frequently-quoted model predictions are 
ft = 0.83 ± 0.08, 0.44 ± 0.03, and 0.25 ± 0.03 for random su- 
perpositions of supernovae in the starburst model, instabilities 
in the accretion disk, and microlensing due to compact lenses in 
the foreground, respectively (Kawaguchi et al. 1998; Hawkins 
2002 ). A tendency of the SF slope to increase with wavelength 
in the observer frame might be indicated in Table|3] (see also 
Schmidt et al. l20T0l their Fig. 3). 

The amplitudes of the SFs shown in Fig.H] generally agree 
with findings in the literature listed in TableHj De Vries et al. 
(2005; their figure 18) present a noise-corrected SF for com- 
bined g and r band as a function of time-lag in the observer 
frame. The last data point in our Fig. [4] corresponds to time- 
lags of ~ 7yr where the de Vries SF predicts log V = -1.08, 
in perfect agreement with log V = -1.0 + 0.1 for g and r in our 
study. De Vries et al. also show (their figure 9) that the time- 
lag shifted SF from Hawkins (2002) agrees with their data. The 
SF derived by Hook et al. ([1994) corresponds to a lower ampli- 
tude of log V = -1.54 at r ~ 7 yr (observer frame). The max- 
imum lag of the data from Bauer et al. ( 12009b is ~ 10 3 d (rest- 
frame), but the amplitudes are not clearly defined at the longest 
time-lags. For log(r r (ii)) = 2.5, Bauer et al. give log V = -1.6 
without specifying the filter band, compared with -1.3 or -1.4 
for g and r in the present paper. However, Bauer et al. empha- 
size that "the exact values for V should not be taken seriously 
. . .because the data are arbitrarily normalized". From the SFs 
presented by Wilhite et al. (120081 see also Rengstorf 2006)) we 
derive log V = (-0.99, -1.17, -1.36, -1.44) for (u,g,r,i) at their 
maximum time-lags log (r r (d)) = 2.7, in good agreement with 
the values (-1.0,-1.15,-1.28,-1.35) found here for the same 
T r . The variability amplitudes derived from the SFs shown by 
Vanden Berk et al. (2004|) are almost the same as in Wilhite et 
al. ((2008). Kawaguchi et al. 119981 ) discuss the SF of only one 
quasar, hence their derived variability amplitude is not necessar- 
ily representative. We emphasize that these results are derived 
from SFs defined in different ways and computed from different 
quasar samples, that cover different ranges of time-lags, and are 
related to different frames (rest or observer). 



H. Meusinger et al.: Spectral variability of quasars in SDSS Stripe 82 



_ 1 1 II 1 1 II 1 1 1 1 1 _ 

— ' ■ ',11 - 


_ l | l l l | l l l | 1 1 1 _ 


_ i | ii i | ii i | ii i _ 


_ 1 | 1 1 1 | 1 1 1 | 1 1 1 . 

- ■ ; - ' i - 


_i i i i i i i i i i i i i_ 

— z — 


~ ■' ' " / '%*• V "'^ "~ 
I 1 I I I 1 I I : 1 i-'.i i 


■ dk^i 




i 1 : i'VI V r- i l-H.V i - 


i ■ ! i i i i 


I 1 I I I 1 i i i 1 . : I 


™...,J 



5 -1 

o 



-22 -24 -26 



-22 -24 -26 
M,- 



-22 -24 -26 
M,- 



-22 -24 -26 



-22 -24 -26 
M,- 



i i | i i i i | i i i i - 



i i i I i i i i I i i i i 



i i i i I i i i i I i i i i 



=* -1 




I i i i I i i i i I : i i 







I I 1 I I : I i-i 1 l 



I I I I | I I I I | I I I I 



- 






_l I I I I I I L 



J_ 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 



J_ 



J I I I I I l_J I I I I I I I I I I I I I I I l_l I L 



12 
redstiift z 



1 2 
redstiift z 



1 2 
redshift z 



1 2 
redshift z 



1 2 
redshift z 



Fig. 5. Variability-absolute magnitude relations (top) and variability-redshift relations (middle) for the five SDSS bands. The bottom 
line compares the measured z gradient of variability (solid) with the z gradient expected from the dependence of V on M;(z) (dashed). 
(Note that the M\ are /f-corrected in the rest-frame, yet the Vk refer to photometric bands defined in the observer frame.) 



4.2. Variability as a function of luminosity and redshift 

The sample-averaged variability measured in a given photomet- 
ric passband is expected to change with redshift because of the 
dependence of variability on absolute magnitude and intrinsic 
wavelength. In addition, there might be an explicit dependence 
of variability on redshift, i.e., V = V[M;(z), A(z); z] and 

dlogV dlogV dM; dlogVd,} dlogV 

— = + + — (3) 

dz 3M; dz 8A dz dz 

where A is the rest-frame wavelength. 

Figure[5] shows the distribution of the variable quasars (x 2 > 
3) in the log V-M\ diagrams (top) and the log V-z diagrams (mid- 
dle), respectively, for the five bands. We overplot the values of 
log(V), where () denotes the median within bins of the width 
0.5 mag for Mj and 0.05 for z, respectively. For luminous quasars 
(here: Mj ^ -25), the average variability decreases with lumi- 
nosity. The relatively small variability of low-luminosity AGNs 
in the z band is expected because of the contribution of the host 
galaxy. The correlation with redshift is less pronounced but a 
slight tendency for larger variability at lower redshifts is appar- 
ent for the u and g band. 

The solid polygon curves in the bottom row of Fig. [5] are 
the gradients d log V/dz (i.e., the left-hand side of Eq.|3]l derived 
from the median relations shown in the middle row. The dashed 



lines represent the fraction corresponding to the M\-z relation 
(i.e., the first term on the right-hand side of Eq.|3]l computed 
from the log(V)-Mj relation from the top line and the (Mj)-z re- 
lation from Fig.[T] In principle, the difference between the two 
curves can be used to derive the wavelength dependence of vari- 
ability provided that there is no significant explicit redshift de- 
pendence of the variability. In particular, that the solid curves 
have generally higher values than the dashed curves at higher 
z indicates that V increase towards lower intrinsic wavelengths. 
However, the V-A relation derived in that way is much too noisy 
to provide useful information. The observed V-z relation is dom- 
inated by the V-M;(z) relation as indicated by the similarity of 
the two curves. A more effective approach to the V-A relation 
will be discussed in Sect. l4.3l below. 

To search for an explicit z dependence, we compare the 
variability measured at different z but the same intrinsic wave- 
lengths for quasars of the same absolute magnitudes. This can 
be achieved with the present sample thanks to the large num- 
ber of quasars and the simultaneous measurements in five bands. 
We select the variable quasars in the narrow absolute magni- 
tude interval M\ = -25 + 0.5, where the coverage in the red- 
shift space is good for z ~ 0.5 . . . 2.2. For each rest-frame wave- 
length A ~ 1000 . . . 4000 A, we computed the redshift intervals 
Zk ± Azk/2, where A is shifted into the u,g,r, and i passbands, re- 
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spectively, i.e., Zk = Ak/A - 1 and Azk = AAk/A where Ak is the 
mean wavelength and AAk the width of the corresponding band 
in the observer frame (k - 1 ... 4). 

The top panel of Fig.[6]shows the mean variability in the red- 
shift intervals around Zk(A) as a function of intrinsic wavelength 
A. At fixed A, the band k corresponds to the redshift Zk(A). At 
wavelengths where there is an overlap of measurements in at 
least two bands, the comparison of the curves provides informa- 
tion on V(z). With a mean value of <dV/dz> = 0.001 ± 0.007 
for 1000 A < A < 3500 A, it is clearly seen (bottom panel) that 
there is no significant explicit variation in V with z for medium- 
luminosity quasars between z ~ 0.5 ... 2. In addition, an increase 
in variability towards shorter wavelengths is clearly apparent in 
all bands (top panel). We note however that the absolute values 
of V(A) are not representative for the whole sample because of 
the V-Mi relation. 



3 




: ooo 



2000 3000 
wavelength \(A) 



Fig. 6. Variability as a function of intrinsic wavelength (rest- 
frame) (top) and gradient dV/dz (bottom) for quasars with 
-24.5 > M\ > -25.5 (same scale for both parts of the diagram). 



4.3. Variability as a function of intrinsic wavelength 



As outlined in Sect. 3.3 variability ratios 2k,i were computed for 



each quasar to measure how variability changes from band to 
band, thus with wavelength. The fundamental advantage is that 
the Qk.\ are expected to be much less (or not at all) contaminated 
by a dependence of variability on luminosity (Sesar et al. 2007). 
The approach is based on the assumption that quasar variability 
is dominated by the same processes at high and low redshifts. 
We argue that this assumption is supported by our finding that 
there is no indication of an explicit dependence of variability on 
redshift in our sample. The variability ratios refer to the same 
intrinsic time-lags of 300 to 600 days for all redshifts. 

In the following, we consider only AGNs with M\ < -21 
and z = 0.2 ... 3 because of the small number per bin out- 
side this redshift interval. As in the previous subsections, we 
remove the relatively small fraction of unusual BAL quasars, 
core-dominated radio-loud quasars, and the few objects with un- 
certain redshifts. The final number of quasars is large enough to 
bin the data into small z intervals of width Az = 0.05. In each of 
the 56 z bins, the sample-averaged ratio (Qk,i(z)) carries informa- 
tion about the variability at five different intrinsic wavelengths 
Ak- The ensemble variability is thus measured at 56 x 5 = 280 



different intrinsic wavelengths. Across the wavelength range 
A ~ 800 . . . 8500A, the formal resolution of the intrinsic variabil- 
ity spectrum is therefore (8500-800)A/280 ~ 28A. The spectral 
resolution can be improved by a factor ~ 2 by choosing narrower 
z bins. 

The results are shown in the third row of Fig. [7] The dots are 
the Qk,i per quasar on a logarithmic scale with k and / referring to 
the SDSS bands (Ak < Ai) shown in the first row of the same col- 
umn. The curves are the log (Qk,i) - z relations where () denotes 
the median in z bins of the widt hs .05. Only quasars with sig- 



nificant variability (x 2 > 3; Sect. 3.1} in both the k band and the 1 
band were included. The number of involved quasars is given in 
each panel. In the high-throughput bands, more than 7000 vari- 
able quasars are available. For most redshifts variability is, on 
average, stronger in the shorter wavelength bands than at longer 
wavelengths. But it is also obvious from Fig.[JJthat the variabil- 
ity ratios display a complicated, nearly undulating dependence 
on z. 

The z-dependence of the variability ratios can be qualita- 
tively understood by comparing with the second row of Fig. [7J 
where the colour indices, mk - mi, are shown as a function of 
redshift. The variability ratios and the colour indices have a 
very similar dependence on redshift, though the scatter is clearly 
larger for the former. This similarity was first noted by Sesar et 
al. (|2007) for the g and r band in the redshift interval 1 .0 ... 1 .6 
and is shown here to also exist for the other bands and across 
a broader z range. The structure in the colour index-redshift di- 
agrams is known to be caused by the strong emission lines. To 
illustrate this effect, we show in the top row of Fig.[7Jthe quasar 
composite spectrum in the observer frame (vertical direction) as 
a function of z (horizontal direction). The spectrum has been 
multiplied with the transmission curves of the corresponding 
SDSS bands. The most prominent emission lines are labelled. 
The colour index Amy decreases when a strong line enters the 
band k and increases at redshifts where the line leaves the band k 
and enters the band /. The effect is most clearly indicated by the 
Mg ii /12800A line where there are no other strong lines around. 
The line dominates the u,g,r,i,z band at z ~ 0.3,0.7, 1.2, 1.7, 2.2, 
respectively, u—g is consequently quite blue at z ~ 0.3 but redder 
at 0.7, g — r is blue at 0.7 but red at 1.2 etc.. 

The variability ratios vary with z in the opposite way to the 
flux ratios. When the flux ratio fk/fi increases because a strong 
line enters the band k, the variability ratio decreases. This can 
be explained only by the assumption that lines are significantly 
less variable than the continuum at the same wavelengths. Strong 
emission lines "dilute" the variability in the corresponding band. 
Hence, the (Qk,i)-z relation contains information on the wave- 
length dependence of the quasar variability and on the relative 
fraction of variability in the continuum and in the lines. In the 
next section, we use a simple parametrized model to quantify 
this relation. There is no hint of a correlation between the de- 
viations of the flux ratios and the variability ratios from their 
median relations. 

The bottom line of Fig. [7J shows that, at a given redshift, the 
variability ratio gk,i also depends on luminosity. Interestingly, 
the gradient d log Qk.\ /d M; varies in a way with z that resembles 
the variation in the variability ratios itself. 



5. Discussion 

5.1. Numerical simulations 

To interpret the results from Fig. [7] quantitatively, we performed 
Monte Carlo simulations of variable quasar spectra. As a first 
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Fig. 7. Top: Spectral coverage of the quasar composite spectrum by two adjacent photometric bands k,l(k = u,g,v,r and I = g,v,r,i,z) 
as a function of redshift z. The next rows (top to bottom) show the corresponding colour indices m\ - m\, the variability ratios (2k,i, 
and the gradients d log Q^\ldM u respectively. Polygons: medians of colour index and log (2k,i, respectively. 



step, the composite spectrum was decomposed into three com- 
ponents: the continuum / CO nt> me emission line portion /ii nes , and 
the host galaxy contamination /hostj^We assumed that the con- 
tinuum is the underlying variable source, the lines are much less 
variable, and the contribution of the host does not vary. 

A simple method was used for the decomposition. The con- 
tinuum / cont was described by the power law fitted to the com- 
posite between Lya and 4000 A, / cont oc .T 048 (Sect.^. The 
power law was extrapolated down to 800 A and a wavelength- 
dependent correction factor cn(A) for the heavy Ly a forest ab- 
sorption was derived shortward of Ly a. The product ch • / CO nt 
was subtracted from the composite and the result shortward of 
4000 A was identified with the UV emission-line portion. The 
result longward of 4000 A is the superimposition of the opti- 



cal emission lines with the host spectrum. The host contribution 
was approximated by fitting by eye a smooth function to the stel- 
lar continuum. After parametrizing the variability as a function 
of wavelength, this rough decomposition enabled us to simulate 
variable spectra and analyse at least the general trends. 



The variability was parametrized by assuming Gaussian- 
distributed magnitude fluctuations, Am, (see Sect. 



3.1 1. The stan- 



The fluxes are always related to wavelength intervals, i.e., f A . 



dard deviation cr Am (A) is described by a power-law dependence 
on the intrinsic wavelength A. Hence, the variable portion of the 
spectrum was represented by f VM (A) = f vm (A) ■ 10"°- 4Am where 
Am = m — m = -2.5 log f vm /f and f vm . was the variable por- 
tion of the composite spectrum. The exponent was computed as 
Am - Am norm cr(/t) where Am norm is a random number for each 
single realization of a variable spectrum. For n realizations of 
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Fig. 8. Comparison of simulations (open symbols) and observations (thin curves and Icr deviations): colour indices (top) and 
variability ratios as functions of redshift for the best-fit model with with small line variability (k - 0.11; second row), the same 
model with large line variability (k — 1 ; third row), and for a model where only the line flux is variable (bottom), respectively. 



the same spectrum, Am n0lm follows a centred and normalized 
Gaussian distribution. 

The aim of the simulation was to compute the variability 
ratios £>k,i as a function of redshift to estimate the parameters 
of the power law as well as the line variability fraction k from 
the comparison of the simulated with the observed data. We 

adopt / var (A) = Ch(A) ■ /contOi) + K ■ fiiaesU) with K = ... 1. 

Hence, a realization of a quasar spectrum is given by f(A) = 
/var( / l) + (1 - k) ■ /ijnesW + /host(^)- The simulations cover the 
same redshift range as the observations. After transforming into 
the observer frame, the "observed" spectrum is multiplied by the 
response functions of the five SDSS filter bands and the total flux 
in each band is transformed into a magnitude. From these data, 
colour indices and variability ratios are computed. 



A parameter study was performed to find the best-fit 
model with 20 realizations of a spectrum per redshift bin. The 
Mersenne Twister MT1993lQ random number generator was 
used, which passes various tests such as the stringent "dieharder" 
test suite for randomness^] The agreement between the simulated 
variability ratios g s ; m and the observed ratios gobs is measured 
by the sum of the (Q s i m - (Gobs)) 2 over all redshift bins, where 
() indicates the average in the redshift bin. We started with the 



3 http://www.dm.unito.it/ cerruti/primi/primigrandi/twister.html 

4 http://www.stat.auckland.ac.nz/dsc- 



2007/abstracts/eddelbuettel06Oct 1 6.pdf 
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assumption that the variability in the lines is negligible (a: = 0). 
The best fit is found for 



/ A \-°- 28 

CTAmW) = 0.06 + 0.07 • r 

V 4000 A/ 



(4) 



The parameter k was subsequently varied; closest agreement was 
found for k ~ 0. 1 . 

The simulated variability ratios of our best-fit model are 
shown in the second row of Fig. [8] The general agreement be- 
tween the observed and the simulated relations is good. Most 
of the features seen in the mean observed relations are closely 
reproduced by the simulations. Exceptions are the low-redshift 
range for the redder bands and the high-redshift range for the 
bluer (z ^ 2 for (2 u ,k an d Z ~ 2.7 for 2g,r)- The former is 
attributed to the rough modelling of the host contribution to 
the spectrum. The latter probably indicates an inconsistency of 
the composite spectrum at shortest wavelengths (A ^ 1100 A) 
with the photometric data producing also a discrepancy for the 
colour indices. It is expected that a more elaborated method for 
the spectrum decomposition and an individual treatment of the 
strongest lines (Sect. 5.2.2[ ) will provide an improved fit. 

The third row of Fig. [SJ shows the results from the same 
model with k = 1. The comparison with model k — 0.1 clearly 
indicates that most of the structure seen in the Qk,i-Z relations 
is due to the small variability of the lines. In addition, the in- 
fluence of the host galaxy is evident in the redder bands at low 
redshifts. For k <: 0.2, the structure in the simulated curves is 
noticeably smeared increasingly with increasing k. Finally, the 
structure reverses for models where only the line component is 
variable (bottom line of Fig. IS). 



5.2. Comparison with Wilhite etal. (2005) 

5.2.1. Colour variability 

On the basis of the analysis of spectra of quasars with repeated 
SDSS spectroscopy, Wilhite et al. (2005) discuss the redshift de- 
pendence of colour differences A(mk - mi) = (m^ - mi)bri g ht - 
(rak - TOi)f a ; nt . There is a change in colour that is a function of 
redshift and the filter used (their fig. 14). As pointed out by those 
authors, the features in these diagrams are related to the relative 
lack of variability in the emission lines. In all colours, quasars 
appear bluer in the bright phase at almost all redshifts. This be- 
haviour is modulated by, with changing redshift, strong emission 
lines being shifted in or out of the passbands. When a strong line 
shifted into the band k is less variable than the continuum, the 
fraction of the line flux to the continuum flux is larger in the faint 
phase than in the bright phase making (m^ - OTi)f a ; nt bluer than 
expected from the power-law continuum and the colour differ- 
ence A(;«k - redder. The features in fig. 14 of Wilhite et al. 
can be directly compared to those seen in the third line of our 
Fig. [7] For example, the first "red bump" in their A(m - g) versus 
z diagram at z ~ 0.35 corresponds to the minimum of Q n g i.e., 
where the variability in u is low compared to g, and the "blue 
dip" at z ~ 0.75 corresponds to a maximum of Q u , g . In this way, 
all dips in the four colour-redshift diagrams in fig. 14 of Wilhite 
et al. can be identified with bumps in the £>k i diagrams in our 
Figs.|7]and|S] 

5.2.2. Emission line variability 



differs significantly from the 10% we find from fitting the vari- 
ability ratios. For k > 0.2, the structures in the simulated Qk.\-z 
relations are remarkably lower and the agreement with the obser- 
vations becomes worse. However, our parameter k is an integral 
quality that is averaged over all lines. 

In this context, we note that even our best-fit model is not 
perfect. In particular, the slight systematic shift towards lower z 
of the structure related to the Mg n line in the simulated curves is 
most likely due to differences in the variability fraction of Mg n 
and the lines in the 3000 A bump. More detailed simulations with 
an individual treatment of the strongest lines will certainly im- 
prove the fit and provide insight into the relative strengths of the 
variability of these lines. 

5.2.3. Continuum variability 

The standard deviation <x in the flux variation as a function of 
wavelength is shown at the top of Fig. [9] for our best-fit model 
from Fig.|SJ The fluxes in the lines were assumed to vary by 
10%, hence the lines are visible in this diagram. For A ^ 2500A, 
the continuum variability is fitted by a power law 



cr(f A ) cc A 



-2 



(5) 



Wilhite et al. (2005) estimate the fraction of variability in the 
lines of from 20% to 30% for Civ, Cm], and Mgn. This result 



in perfect agreement with the slope of the geometric mean dif- 
ference spectrum from Wilhite et al. ([2005). This means in par- 
ticular that the variability spectrum is steeper (i.e., bluer) than 
the flux spectrum. 

Wilhite et al. ([2005 ' found that the ratio of their arithmetic- 
mean composite difference-spectrum to the arithmetic-mean 
composite quasar-spectrum is almost one for A S 2500 A and 
the ratio strongly increases towards shorter wavelengths at A 55 
2500 A. The ratio of the standard deviation in our best-fit model 
to our arithmetic-mean composite is shown in the bottom panel 
of Fig. [9] We confirm that this ratio is approximately constant in 
the wavelength range between the the Mg n line and A ~ 4000 A. 
However, there is a dip shortward of Mg n in both our Fig. ^(bot- 
tom) and the corresponding fig. 13 of Wilhite et al., which leads 
to the impression of an abrupt increase shortward of ~ 2400 A. 
In our model, this dip is attributed to the very broad conspicuous 
"contamination" feature from ~ 2300 to 4000 A (see Vanden 
Berk et al. l200TT) . This so-called small blue bump (3000 A bump) 
consists of a large number of weak Fe n emission lines blend to- 
gether to form a pseudo-continuum above the intrinsically emit- 
ted continuum, superimposed on the Balmer continuum emis- 
sion (Wills et al. [19851 Vestergaard & Wilkes 1200Tb . We argue 
therefore that the variability increases rather continuously to- 
wards short wavelengths for A 4000 A. At longer wavelengths, 
our ratio decreases dramatically, in contrast to that of Wilhite et 
al.. However, this difference is expected since our sample is dom- 
inated by low-luminosity AGNs at low redshifts and the contri- 
bution from the host galaxies [A <: 4000 A) to our composite 
spectrum is thus much larger than in Wilhite et al.. 

5.3. Luminosity dependence of variability ratios 

The behaviour of the luminosity dependence of the variabil- 
ity ratios as a function of redshift (bottom row of Fig. [7) can 
be qualitatively understood as follows. We assume that F C om,k 
and Fii ne ,k are the total fluxes of the continuum and the emis- 
sion lines, respectively, in the band k. If the lines are far less 
variable than the continuum, the fraction of the continuum flux 
Pk = F contX /(F conUk + Fi ine ,k) is a measure of the strength of 
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Fig. 9. Top: Wavelength dependence of the standard deviation 
cr of the flux density fx for the best-fit model with k — 0.1. The 
dashed line corresponds to the power law with slope -2. Bottom: 
Ratio cr(f /t )/f A where /, is the flux density of the composite 
spectrum normalized at 4000A. 



variability at a given redshift. The variability ratio Qk.i is then 
expected to scale with the ratio p^/pi- For an average quasar 
spectrum and at given z, the luminosity scales with the total flux 
^tot = -Fcont + -Fiine in a given band. Hence 



G = 



dlog g kj] F tot ,k dgk,i 



dMi 



2k,l dF toti k 



For simplicity, we consider the case where there is only one 
strong emission line (e.g., Mgn) in band k for z ~ Zi and in 
band 1 for z~ Zi > Z\ . Hence 



/ r con t \ j / . d/* line \ c 

Pk = [-=— . A = l, G k oc-g M+ — forz-z, 

pi _ / cont\ pk = l j Gi oc -gk,l (l + -TTT^-) forz~Z 2 - 



The gradient G depends directly on the variability ratio and, 
in addition, on the relationship between the line flux and the 
continuum flux in the band. The equivalent widths of the broad 
UV emission lines is known for a long time (Baldwin 119771 
Osmer et al. 119941 ) to be inversely correlated to the luminosity of 
the underlying continuum in a nonlinear way. This well-known 
Baldwin effect appears in two flavours: (i) a global relationship 



that describes fluctuations from object to object and (ii) an in- 
trinsic relationship that describes how the line flux changes when 
the continuum of the same object changes because of variability 
(Kinney et al. 1990; Pogge & Peterson 1992). In addition, how- 
ever, the variation in G with z also reflects the response functions 
of the SDSS bands because the ratio of the line flux to the con- 
tinuum flux in the band depends on the throughput at the wave- 
length of the line. To summarize, the G-z relation (bottom line 
of Fig. [7]) reflects the Qy\-z relation, which is, however, modified 
by both the Baldwin effect and the filter response curves. 

5.4. Special quasar types 

Among the 8744 objects in our final quasar sample, there 
are 54 quasars with unusual spectra and 454 quasars were 
identified with FIRST radio sources. Adopting the criterion 
Ri s log (^radio/^i) > 1 (Sect.|2]i, we identified 381 radio-loud 



quasars in our sample. In Fig. 10 the optical colours and vari 



ability properties of both subsamples are compared with those 
of our standard quasar sample. As was noted by Richards et al. 
(2001) and Ivezic et al. (2002), the median colours are slightly 
redder and the fraction of red outliers is higher for the radio-loud 
subsample than for the standard sample, where the most extreme 
outliers are seen in the bluer passbands. A substantial fraction of 
the radio-loud subsample are also outliers in the Qk,\-z diagrams. 
In particular, there is more than a dozen quasars with low 2k,h 
i.e., Vi > 2V t . 

For the subsample of quasars with unusual spectra, the de- 
viation from the standard sample is much stronger. The colour 
indices u - g and g - r are much redder than for both the stan- 
dard sample and the radio-loud subsample. These differences are 
again most pronounced for the bluer passbands and vanish for 
i - z. This result is not surprising. Both dust absorption and ex- 
tended and/or overlapping BAL troughs in the UV, which are two 
defining features of this classification, have the effect of redden- 
ing especially the colour indices of the blue passbands. In this 
context, the scatter in the variability ratios might be expected to 
be explained by the diversity of the fraction of the emission-line 
flux in strong and unusual BAL quasars. 

The fraction of radio detections among the unusual quasars 
is 17% (9/54) and the fraction of radio-loud quasars is 
11% (6/54). The corresponding values for the whole sam- 
ple are 5.2% (454/8744) and 4.4% (381/8744). This dif- 
ference is, of course, due to SDSS having targeted FIRST 
sources for spectroscopy even when they were not selected 
as quasar candidates because of their optical colours. The 
three unusual quasars with the highest radio-loudness param- 
eters are SDSS J00005 1.56+001202.5 (Ri = 1.89), a high- 
redshift (z = 3.859) BAL quasar with broad Lyo-, Siiv, and 
Civ absorption troughs, SDSS J024224.02+010452.5 (R, = 
1.52), a medium-redshift (z = 2.431) quasar with appar- 
ently weak emission lines and associated absorption lines, and 
SDSS J022337.76+003230.6 {R, = 1.77), a weak-line quasar 
at z ~ 3.05. The first two are strongly variable in the optical, 
whereas SDSS J022337.76+003230.6 exhibits a remarkably low 
level of variability in the u and g bands. 

As a main result of Sects. [BTT] and |5.2[ we have seen 
that the flux in the emission lines is considerably less vari- 
able than the underlying continuum. A relatively low vari- 



ability V g is thus expected at z 



3 where the strong 



Lya+Nv lines fall into the g band (see also Fig.|5](. However, 
the value of V u for SDSS J022337.76+003230.6 is a factor 
of 20 smaller than the mean value at z = 3 and, sur- 
prisingly, SDSSJ022337.76+003230.6 is a member of the 
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Fig. 10. Colour indices, - m\, (top) and variability ratios, <2k,i, (second row) as a function of redshift and variability index, Vk, 
as a function of absolute magnitude, M\, (third row) for radio-loud quasars (dots) and quasars with unusual spectra (squares and 
arrows), respectively. Curves and hatched areas as in Fig. [8] filled squares in the third row: unusual quasars with weak emission 



lines (Fig. 1 1 left panel). Bottom: Normalized distributions of log Vk for the unusual quasars (open squares and dashed lines), the 



radio-loud quasars (dots and dotted lines), and the standard sample (solid). 



rare class of weak emission-line quasars (WLQs; Diamond- 
Stanic et al. 2009). WLQs are defined as quasars with ex- 
tremely weak or undetectable emission lines of rest-frame 
equivalent width EW(Lyar + Nv) < 10 A compared to 
~ 50... 100 A for typical SDSS quasars (e.g., Fan et al. 
119991 Diamond-Stanic et al. 2009 and references therein). For 
SDSS J022337.76+003230.6, Diamond-Stanic et al. measured 
EW(Lyo' + Nv) = 8.6 A. Quasars with such weak lines might 
be expected to exhibit larger variability in the passbands of 



the lines than typical quasars. This is clearly not the case for 
SDSS J022337.76+003230.6. Altogether the Diamond-Stanic 
WLQ sample contains four quasars in S82, three of them being 
members of our quasar sample where two have small variability 
(SDSS J022337.76+003230.6 and SDSS J00542 1 .42-010921 .6) 
and the third (SDSS J025646. 56+003858. 3) has values close to 
the medians. 

If the weak lines are due to the Baldwin effect, WLQs might 
be expected to populate the high-luminosity end of the quasar 
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luminosity function, which is obviously not the case (Shemmer 
et al. 120091 and references therein). The first interpretation of the 
nature of WLQs involved continuum boosting (Fan et al. 1999) 
as in local BL Lac objects. However, the rest-frame 0.1-5/im 
spectral energy distributions of WLQs does not differ from those 
of normal quasars and variability, polarization, and radio prop- 
erties differ from those of BL Lacs (Diamond-Stanic et al.2009 ). 
In addition, Shemmer et al. (2009) pointed out that it would be 
difficult to explain the lack of a large parent population of X- 
ray and radio bright weak-lined sources at high z if WLQs were 
the long-sought high-redshift BL Lacs. Shemmer et al. ([2006) 
argued against the possibility that WLQs are quasars with con- 
tinua amplified by microlensing. This interpretation is also ques- 
tioned by the low variability found in the present study based on 
intrinsic timescales of 1 to 2 yr, which are typical timescales for 
quasar microlensing. 

Among the quasars classified here as unusual, no additional 
one matches the classical criterion of WLQs (i.e., EW(Lya + 
Nv) < 10 A). We note that our quasar sample is biased 
against spectra that have only a featureless continuum (Sect.|2|). 
However, we selected (by eye) another five quasars for which 
the emission lines are substantially less pronounced than in 
the composite spectrum but where there is no clear-cut evi- 
dence of BAL troughs. Their spectra are shown in the left panel 
of Fig.fTT] along with the WLQ SDSS J022337.76+003230.6. 
SDSS J220445.27 +003141.8 (z = 1.353) is one of the two "mys- 
terious objects" from Hall et al. (2002). The interpretation of 
SDSS J010605.07-001328.9 is unclear but the spectrum resem- 
bles the unusual F eLoBA L quasar VPMS J134246.25+284027.5 
(Meusinger et al. 12005). In the log V-M\ plots in Fig. 10 these 



quasars are marked as filled squares. All of them vary by less 
than the general median dispersion. 

Apart from the puzzling WLQs, the variabilities of the un- 
usual subsample are systematically smaller than those of the 
standard sample. The bottom row of Fig. 10 displays the dis- 
tributions of log Vk for the unusual quasars, the radio-loud 
quasars, and the standard sample in the absolute magnitude 
range -27 < M{ < -23. There are no substantial differences 
between the radio-loud subsample and the standard sample. 
For the unusual quasars, however, we find that 24.5% (13/53) 
of the quasars with M, < -23 have V < 0.01 in each of 
the bands u,g,r, and i, compared with 6.6% (483/7334) for 
the standard sample in the same M; range. Among these 13 
low-variability unusual quasars, there are in particular all six 



WLQ-like objects from the left panel of Fig. 1 1 The spectra 
of the remaining 7 quasars are shown in the right panel of 
Fig. 1 1 most of which have strong and/or unusual BAL struc- 
tures. The properties of both SDSS J025858. 17-002827.1 and 
SDSS J024543. 36+000552.2 are most likely caused by many 
low ionisation Fen absorption troughs (FeLoBALs). It is re- 
markable that in all these spectra there is in particular a the lack 
of strong emission lines. 

Why is the variability of the unusual quasars smaller than 
the median values, even though a lack of broad-line emis- 
sion is expected to result in stronger variability? Wilhite et al. 
(2008) proposed that quasar variability is inversely related to 
the Eddington ratio Lboi/^Edd, and Shemmer et al. Q.Q09 ) sug- 
gested that WLQs might be quasars with high Eddington ratios. 
Ai et al. (2010) found by comparing broad-line Seyfert 1 and 
narrow-line Seyfert 1 type quasars that the variability amplitude 
correlates with the relative strength of emission lines and anti- 
correlates with the Eddington ratio. If Lboi/^Edd is an ( at least 
crude) estimate of the accretion rate, both weak emission lines 
(but strong Fe n pseudo-continuum) and small variability ampli- 
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Fig. 11. Left: Unusual quasar spectra with weak emission lines. 
Right: Another 7 unusual quasars with weak variability {right). 
In each panel, the composite spectrum is overplotted as smooth 
curve (arbitrarily scaled). 



tudes are the consequence of strong accretion. An anticorrelation 
with variability is plausible because, if the accretion rate is high, 
the region emitting the continuum observed in a given pass- 
band is larger and the variability amplitudes are consequently 
smaller than in the case of a low accretion rate (Ai et al. 120101 ). 
Wilhite et al. (2008 ) proposed that the accretion rate is high and 
the variability is small for optically young quasars, i.e., when 
they become visible in the optical. On the other hand, FeLoBAL 
quasars have been suggested to represent optically young AGNs 
expelling a thick shroud of gas and dust from the cocoon phase 
(e.g., Voit et al. 119931 Farrah et al. 120101 ). In this context, the 



higher fraction of low-variability quasars among the unusual ob- 
jects might be an indication that the unusual subsample contains 
a higher fraction of young quasars. However, the present sample 
is quite small and inhomogeneous. A more detailed investiga- 
tion of the properties of quasars with atypical spectra in S82 is 
necessary. 

6. Conclusions 

We have compiled a sample of nearly 9000 quasars from the 
SDSS Stripe 82 with redshifts z ~ 0.1 ... 5 and absolute magni- 
tudes M; 20 ... - 30 that have multi-epoch photometry in the 

LMCC (Bramich et al. (2008 ). The fraction of variable quasars 
in the sample is found to be 93%, 97%, 93%, 87%, and 37% 
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for the u,g,r,i, and z band, respectively. The SDSS spectra of 
all these objects were studied individually. For the wavelengths 
range from Lya to ~ 4000 A, our arithmetic -mean compos- 
ite spectrum is similar to the composite spectrum from Vanden 
Berk et al. (2001). At longer wavelengths, our composite is more 
strongly contaminated by stellar radiation from the host galaxies. 
The continuum is well fitted by a power law with a v = -0.48. 

After removal of strong photometric outliers caused by "bad 
epochs" in the LMCC lightcurves, first-order structure functions 
(SFs) were computed and used to analyse the ensemble variabil- 
ity of the "normal" quasars (excluding unusual BAL quasars and 
the known core-dominated radio-loud quasars). The following 
main conclusions are drawn: 

1 . The SF as a function of the time-lag can be fitted by a power 
law with the index j3 ~ 0.31 + 0.03 for rest-frame lags 
r = 100 . . . 1000 d. Our data confirm the existence of anti- 
correlations of variability with both luminosity and redshift 
Z. The trend with z can be interpreted as a consequence of the 
dependence on both luminosity and intrinsic wavelengths. 
There is no indication of an additional redshift dependence. 

2. The dependence of variability on rest-frame wavelengths can 
be established by studying the ratio of the variability in two 
different passbands as predicted by Sesar et al. (2007). The 
redshift dependence of the variability ratios, £>k,i, resembles 
that of the corresponding colour indices. The Q^yz relations 
indicate that variability is almost always stronger in the blue 
passband than in the redder one. However, the relations are 
modulated by the appearance of strong emission lines in ei- 
ther one band or the other. 

3. By means of simple Monte Carlo simulations of spectral 
variability, we have shown that the observed Q^\-z rela- 
tions can be described by assuming a variability spectrum 
f(//t) K ^ 2 for me continuum (<x: standard deviation) and 
a low percentage of ~ 10% of variability in the emission 
lines. Quasars are mostly bluer in their bright phase but this 
depends on both z and the passbands. The low variability in 
the lines is related to the intrinsic Baldwin effect. These re- 
sults, based upon the photometry of ~ 8000 quasars, confirm 
the findings of Wilhite et al. (2005) for 315 quasars with re- 
peated SDSS spectroscopy. 

4. The microlensing scenario for quasar variability is not ex- 
cluded by the observed slope of the SF. On the other hand, 
the derived fluctuation spectrum is consistent with the as- 
sumption that the dominant fraction of the optical/UV quasar 
variability comes from the accretion disk as shown by 
Pereyra et al. ( 120061 ). It remains to be seen whether the mi- 
crolensing scenario predicts the observed ratio of variability 
in line flux to continuum flux. But prior to that, more detailed 
information on the variability of the various line components 
needs to be extracted from the data. 

5. For quasars with spectra that deviate significantly from the 
quasar composite spectrum, not only the colour- z but also the 
Gk.rz relations differ from those of typical quasars. Unusual 
quasars with weak or undetectable emission lines tend to ex- 
hibit smaller fluctuations than normal quasars. This trend, 
obviously opposite to the "dilution" of the observed vari- 
ability by the emission line flux, might be an indication that 
a significant fraction of these unusual quasars accrete with 
high Eddington ratios. 
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